IR divergence does not affect the gauge-invariant curvature perturbation 
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We address the infrared(IR) divergence problem during inflation that appears in the loop corrections to the 
primordial perturbations. In our previous paper, we claimed that, at least in single field models, the IR divergence 
is originating from the gauge artifact. Namely, diverging IR corrections should not appear in genuine gauge- 
invariant observables. We propose here one simple but explicit example of such gauge-invariant quantities. 
Then, we explicitly calculate such a quantity to find that the IR divergence is absent at the leading order in the 
slow-roll approximation for the usual scale invariant vacuum state. At the same time we notice that there is a 
subtle issue on the gauge-invariance in how to specify the initial vacuum state. 
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I. INTRODUCTION 

The precise measurements of the primordial fluctuation 
provide us with valuable information of the early universe. 
It is widely accepted that the primordial fluctuation originates 
from the quantum fluctuation of the inflaton field. In the last 
decades, it has been recognized that the perturbation theory 
in the inflationary universe might break down because of the 
infrared(IR) divergence from loop corrections [1]. (For a re- 
cent review, see Ref. [2].) During inflation, massless fields are 
known to yield the scale invariant spectrum P(k) ex 1/fc 3 at 
linear order. These fields contribute to the one-loop diagram 
with the four point interaction as J d 3 fc/fc 3 , which leads to the 
logarithmic divergence. 

It has been an issue of debate whether the IR divergences 
are physical or not. If it were really physical, there might be 
a possibility that loop corrections are observable. The pos- 
sibility of secular growth of IR contributions has also been 
studied motivated as a possible explanation of the smallness 
of the cosmological constant [3], However, before we discuss 
implications of the IR effects, we need to carefully examine 
whether the reported IR divergences are not due to careless 
treatment. (See also Ref. [4]) 

In our previous work [5], we pointed out the presence of 
gauge degrees of freedom in the frequently used gauges such 
as the comoving gauge and the flat gauge, and these gauge 
degrees of freedom are responsible for the IR divergences. 
There, we have shown that, if we fix the residual gauge de- 
grees of freedom, the IR divergences automatically disappear 
in the single field model. (Multi-field case was discussed in a 
separate paper [6].) 

In this paper, we reexamine our previous argument. If the 
IR divergences are really due to the residual gauge degrees of 
freedom, those divergences should disappear if we evaluate 
genuine gauge-invariant quantities even though we do not fix 
the residual gauge. Hence, the demonstration of the genuine 
gauge-invariance of computed results would be necessary to 
obtain reliable predictions, which are to be compared with ob- 
servations. In this brief report we demonstrate more explicitly 
that the IR divergence in the curvature perturbation can be re- 
moved by looking at such genuine gauge-invariant quantities. 
In doing so, we shall notice that the choice of initial vacuum 
state is restricted. In the slow roll limit, we find that a natural 



vacuum state is surprisingly limited to the Bunch-Davies vac- 
uum state as far as we consider Gaussian vacuum state for the 
inflaton perturbation in the flat slicing at the initial time. 



II. BASIC EQUATIONS 

We consider the single field inflation model whose action 
takes the form 



S 



-y 1 / V=9 [R - S^<MV - 2V((b)]d 4 x , (1) 



where M p \ is the Planck mass and the scalar field $ was 
rescaled so as to be non-dimensional as in Ref. [5]. Substi- 
tuting the metric form 

ds 2 = ~N 2 dt 2 + h ij (dx i + N l dt)(dx j + N j dt) , (2) 

the action becomes 



S = 



M, 



pi 



Vh 



N S R- 2NV{(f)) + —(EijE 11 - E 2 ) 



+ iV (0 ~ N%di(/) > ~ KV'fyWrf 



where Eij and E are defined by 



d 4 a; , (3) 



E^-^-D^-D^), E = h- Eij , (4) 

In this paper we work in the comoving gauge, defined by 
p = 0. We decompose the spatial metric as 



e a(/H-C) 7tf = ^0 [ e *r] 



(5) 



where a :— e p is the background scale factor, and tr[<$7] = 1, 
Using the degrees of freedom of choosing the spatial coordi- 
nates, we further impose the gauge conditions d l Sjij = 0. 

Varying the action with respect to N and N l , we obtain the 
Hamiltonian and momentum constraints as 

S R -2V- N- 2 (E ij E i:j - E 2 ) - N- 2 <p 2 = , (6) 
D 3 [N- 1 (E^-P l E)]=Q. (7) 



For convenience, we factorize the scale factor from the metric 

as 

ds 2 = e 2p [-(iV 2 - N t N l )dr/ 2 + 2N l dr ) dx l + /jydzW], 

where we have defined hij := e~ 2p hij = e 2 ^7y, Ni = e p N, 
and N l := h^ Ni = e p N t . Expanding the perturbations, 
Q = SN(:= N - I), N,(, and <5 7ij as Q = Qi + ±Q 2 + 
• • • , the zeroth-order Hamiltonian constraint equation yields 
the background Friedmann equation: 



6p' A = <t>' z + 2V{<t>) 



(8) 



where a prime " ' " denotes the differentiation with respect to 
the conformal time 77, and V(<j>) :— e 2p V((f>). The constraint 
equations at the linear order are obtained as 



irSNi - 3p'([ + 9 2 Ci + p'd l N lA = , 



(9) 



49, (p'SN, - CO - d 2 N iA + diVftj,! = . (10) 
The higher-order constraints can be obtained similarly. 



III. ORIGIN OF IR DIVERGENCE 

A. Influence of boundary conditions 

Solving the constraints (6) and (7), we can express SN and 
N in terms of the curvature perturbation (. Here we stress 
that the constraints (9) and (10) are elliptic-type equations, 
which require boundary conditions to solve. Below we study 
the relation between the residual gauge degrees of freedom 
and the integration constants appearing in the general solution 
of Eqs. (6) and (7). For illustrative purpose, we consider the 
linear order, in which constraints are given by Eqs. (9) and 
(10), but the extension to higher orders proceeds in a similar 
manner. 

Eliminating SN\ from Eqs. (9) and (10), we obtain 



i-^-l di&tij,! - d 2 N lA 



V 



i,'2 



+ 2^-d^ - %d 4 d 2 Ci = . (ii) 
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V 



Taking the divergence of Eq. ( 1 1 ), the longitudinal part of A^i 
is solved as 

0*^,1 (3) = I^CiO*) - VCiW + Fx(x) , (12) 
Ip z p 

where F\ (x) is an arbitrary solution of the Laplace equation. 
Using Eqs. (11) and (12), N^ is integrated to give 

*iAx) = di(J^d-%(x)-±c 1 (x) 

1 - ^-) d l d- 2 F l (x) + Gi(x) , (13) 



where Gi(x) is an arbitrary vector that satisfies d 2 Gi(x) = 0. 
Here we assumed that the inverse Laplacian d~ 2 is defined 
uniquely. Comparing the divergence of Eq. (13) with Eq. (12), 
we find F x {x) = {V/Ap' 2 )d l Gi{x). Substituting Eq. (12) 
into Eq. (10), the lapse function is easily obtained. To con- 
clude, the lapse function and the shift vector are given by 






2p> 



?' 



(14) 



Introduction of a vector Gi(x) also modifies the evolution 
equations of the curvature perturbation ( and the transverse- 
traceless perturbation Sjij. Substituting Eqs. (14) and (15) 
into Einstein equations, we obtain 

(Ci - A)" + 2(lnz)'(d - A)' - a 2 (Ci - A) - , (16) 

where we have defined z :— e p (f>'/p' and f\{x) := 
2 J drjd z Gi(x), which implies d 2 f\(x) — 0. While the evo- 
lution equation of 5-fij is obtained as 

^-\ didj(d- 2 d 2 -i)Ci + ^did^Gu 

-e- 2p d 11 {e 2p {d l G :j +d 3 G l 

- ±d z d 3 d- 2 d k G k - \d k G k 5 l3 )} = . (17) 

This equation is consistent with the transverse traceless condi- 
tion of Sjij. Namely, the trace and divergence of Eq. (17) are 
automatically satisfied. These evolution equations of £i and 
Sjij y i reduce to the ordinary well-known linear perturbation 
equations by simply setting Gi = 0. 

When we consider the universe with infinite volume, the 
solution of SNi and JV^i would be specified uniquely by the 
requirement that all quantities are regular at the spatial infin- 
ity. In this case, we have not room to introduce Gi (x). In con- 
trast, when we concentrate on a finite region of the universe, a 
variety of solutions are allowed especially due to the presence 
of homogeneous solutions of Laplace equation. In this case, 
as presented by the first term in the second line of Eq. (17), 
it also happens that the evolution equation for the transverse 
traceless mode d-jij is contaminated by the longitudinal mode 
C, even at the linear order. 

B. Gauge degrees of freedom 

The ambiguity in the choice of the vector Gi indicates the 
presence of residual gauge degrees of freedom. Here, we 
show explicitly that Gi(x) represents the residual gauge de- 
grees of freedom that remain even after we specify the gauge 
by 8(j> = and the conditions (5). Since the condition 6<f) = 
completely fixes the temporal gauge, we are only allowed to 



change the spatial coordinates. Under the change of the spa- 
tial coordinates x % — ► x 1 + 8x l , the metric functions transform 



Ni,i(x) = Ni,i(x) - dx'i , 



(18) 
(19) 
<J7»j,i = (fyij.i - (9jfej + Sj&e* - ld k Sx k S lj ) , (20) 

where 6xi := SijSx^ . We associate a tilde with the perturbed 
variables when we set Gj ^ to discriminate them from those 
in the gauge with G, = 0. Recalling that the lapse function 
remains unchanged, the comparison between Eqs. (14) and 
(19) gives 

d l 5x[ = -(3/4)d t G l (x) . (21) 

Imposing the transverse condition of Sjij in Eq. (20) yield 

d 2 5xi = -(l/3)d i d j 6x j . (22) 

Comparing Eq. (15) with Eq. (18) together with Eqs. (21) and 
(22), Sxi is solved as (See Ref. [7].) 



Sxi(x) 



dr/G,(x) + i fdndid-^Gj^x) 
+ \J^T J^d l &G 3 (x)+H l (x) 

+ J±L d i Hi{xh 

where we introduced a vector Hi(x) that satisfies 
?>d 2 H l {x) + d l d J H 3 {x) = {). 



(23) 



(24) 



Substituting Eq. (23) into Eqs. (19) and (20), we find that spa- 
tial components of metric perturbation transform as 



Ci = Ci + lldnira.-Urii;. 



(25) 



*fo,i = *7tf,i + dV {2%G j) - \{didjd- 2 + 5 l3 )d k G k ) 



- l -J'^J&r 1 d i d j d k G k - 2 {d {i H j} - \&H k 8y 



% 2 d (i H. 



pf 



(i n 0)- 



(26) 



To remedy the IR divergence due to the gauge modes, in 
our previous work [5, 6], we proposed to use the local gauge 
condition such that the effects from the causally disconnected 
region are shut off. The local gauge condition can be achieved 
by choosing the function G{ (x) appropriately. To remove the 
IR divergence explicitly, the residual gauge must be fixed so 
that £ measures the deviation from the average value taken 
over a certain local region 0. In Ref. [5], we have shown that 
the IR corrections become finite by adopting this local gauge. 
In this gauge IR contributions are effectively cut off at the 
length scale of the size of 0. However, one can change the 
choice of 0, which leads to a finite shift in the gauge function 
Gi(x). As a result, the resulting correlation functions for ( 



depend on the artificial choice of 0. This prescription will be 
sufficient just to show that the loop corrections are finite in 
some gauge, but it will not be sufficient to evaluate the size of 
loop corrections. In the present work, we give an example of 
genuine gauge-invariant quantities that do not depend on how 
we fix the residual gauge degrees of freedom. The expected 
correlations in the CMB temperature fluctuations should be 
also such gauge-invariant quantities. 



IV. GAUGE INVARIANCE AND IR REGULARITY 

A. Construction of gauge-Invariant variables 

In this subsection, we construct an example of genuine 
gauge-invariant quantities. Since the temporal slicing is al- 
ready fixed, it is sufficient to address the gauge-invariance 
under the residual gauge transformation of the spatial coor- 
dinates (23). Genuine gauge-invariance is equivalent to com- 
plete gauge fixing. In this sense, if we give appropriate bound- 
ary conditions for the lapse and shift, the residual gauge de- 
grees of freedom can be completely fixed. However, here 
arises a difficulty in fixing the gauge, because we need to 
remove all arbitrariness regarding the choice of coordinates, 
such as the choice of the local region 0. 

When we consider the universe with infinite volume, the 
conventional gauge-invariant perturbation theory is formu- 
lated using particular combinations of perturbed variables in- 
variant under an arbitrary gauge transformation. The con- 
struction of such gauge-invariant variables is possible thanks 
to the absence of the ambiguity in the inverse Laplacian un- 
der the assumption that the perturbations remain finite at the 
spatial infinity. However, the same procedure does not work 
when we try to use only the information contained in the lim- 
ited area 0. Under this restriction, genuine gauge-invariant 
quantities cannot be constructed by the combination of local 
quantities. 

Keeping these in mind, we consider n-point functions of the 
scalar curvature S R. The scalar curvature S R does not remain 
invariant but transforms as a scalar quantity under the change 
of spatial coordinates. If we could specify its n arguments 
in a coordinate-independent manner, the n-point functions of 
S R would be gauge-invariant. This can be partly achieved by 
specifying the n spatial points by the geodesic distances and 
the directional cosines from one reference point. Although we 
cannot specify the reference point in a coordinate independent 
manner, this gauge dependence would not matter as long as 
we are interested in the correlation functions for a quantum 
state that respects the spatial homogeneity and isotropy of the 
universe. 

As an example, we consider the two point function of S R, 
whose arguments X A {A = 1, 2) are specified by solving the 
three-dimensional geodesic equation from A = to 1 with 
the initial "velocity" &X A {X) / &\\\ = q = x a starting with the 
origin. We denote X A (\ = 1) simply by X A , which can be 
perturbatively expanded as X A := x\ + 5x A . x A can be 
understood as the background value of the points. Then, the 



expectation value of a product of 
W( V , x\) := 'R( v , x\ + 5x\) 



00 A *! A ^ n 

z — ' n\ 

n=0 



should be gauge-invariant. 



B. IR regularity of genuine gauge-invariants 



some manipulations, the possibly divergent terms in ( 9 R 9 R)4 
can be summed up as 



(W(x) g R(y)U <x (£)J d(\og k)[D 2 u k (x)u* k (y) 

+2Du k (x)Dul(y) + u k {x)D 2 u* k (y)} 
+ (cc), (31) 



The genuine gauge-invariant quantities we introduced 
should be finite in the local gauge, since the field itself is con- 
structed to be free from IR divergence in this gauge by con- 
struction. However, since these quantities are really gauge- 
invariant, they should be finite even if we send the size of 
to infinity. This limit is supposed to agree with the case when 
we calculate them in the infinite volume [7]. In the rest of 
this paper we study the regularity of the gauge-invariant two 
point function introduced above in the conventional global 
gauge [8], focusing on the loop of longitudinal modes at the 
lowest order in the slow-roll approximation. The spatial cur- 
vature S R is given by 



S R 



-2e- 2(p+ ' :) (2d 2 C + d l Cd l O 



(28) 



We expand the gauge-invariant spatial curvature as 9 R = 
9 R\ + 9 i?2 + • • • , in terms of the interaction picture field oper- 
ator. Here the subscript denotes the number of the operators. 
Then, one-loop contribution to the two point function starts 
with the quartic order, 



{ 9 R 9 R)i ■- ( ff i?i ff i? 3 > + ( 9 R2 9 R2) + ( 9 R 3 9 Ri) 



(29) 



Using the fact that the IR divergence arises only from the con- 
traction between interaction picture fields without any deriva- 
tive, we keep only the terms that are possibly divergent. For 
instance, the terms that include more than two interaction pic- 
ture fields with spatial or temporal derivatives do not yield 
divergences. 

In order to obtain the expression for S R, we use the fact that 
in the flat slicing the interaction Hamiltonian is totally sup- 
pressed by the slow roll parameter [8]. We therefore solve 
the non-linear evolution of perturbation in the flat gauge, and 
transform the results into the 6<f> = gauge. Then, the contri- 
butions at the lowest order in slow-roll expansion to ( arises 
only from the non-trivial gauge transformation given by [8] 



C — Cn + 2 P '^ n ^ n + 2p' 2 



(30) 



where £„ := —p'ipi/(f>' and the abbreviated terms are higher 
order or irrelevant for the divergences. Here ipj is the interac- 
tion picture field corresponding to the inflaton perturbation in 
the flat gauge. 

Solving the spatial geodesic equation on a 77 =constant hy- 
persurface, we obtain X 1 = e~^x l + • • • , where the terms 
with spatial differentiations were abbreviated since they do 
not contribute to the possible IR divergence in ( 9 R 9 R)4. After 



where u k (x) is the positive frequency function of the field £„ 
multiplied by fc 7 / 2 , and D is an operator defined by D := 
rjdrj + x' l di + 2. Since (£ 2 ) is IR divergent, the regularity is 
maintained only when the integral in (31) vanishes exactly. If 
we use the scale-invariant property of mode functions of the 
Bunch-Davies vacuum, the derivative operator D acting on 
life is replaced with d\ og k ■ Then, this integral becomes total 
derivative and vanishes. The usage of the invariant distance is 
requested also in Ref. [9], where the SN formalism is utilized. 
It is intriguing that, while the method is different, they have 
arrived at the same conclusion as we have. 

One may think it strange to require the scale-invariance, be- 
cause genuine gauge-invariant correlation functions should be 
finite independent of the choice of initial vacuum state. This 
finiteness is guaranteed as long as we compute it in the local 
gauge. This gauge dependence stems from that of the initial 
vacuum state. In Ref. [5] we implicitly specified the initial 
vacuum state so that the interaction picture field is identical to 
the Heisenberg field on the initial surface. Namely, the inter- 
action is assumed to be turned off before the initial time. Since 
the effect of the coordinate transformation on the perturbation 
variables is non-linear, this identification between the Heisen- 
berg field and the interaction picture field has different mean- 
ing depending on the choice of gauge. Therefore, in order to 
obtain genuine gauge-invariant results, we need the principle 
to specify the vacuum state in a gauge-invariant manner. If we 
specify the initial state so that the interaction picture field is 
identical to the Heisenberg field at the initial time in the local 
gauge, the selected quantum state depends on the choice of the 
local volume 0. This situation is unsatisfactory. Moreover, 
such a choice of vacuum will violate the three-dimensional 
translational invariance. Hence, it is much favoured to spec- 
ify the initial vacuum state in the global gauge. As we have 
seen above, at the lowest order of slow roll approximation, 
the scale-invariant Bunch-Davies vacuum state is such a vac- 
uum state that is free from IR divergences. Since the origin 
of the possible IR divergence is confined to the initial vacuum 
state, we may say that, if genuine gauge-invariant correlation 
functions of fluctuations are free from divergence at the initial 
time, it remains so in the course of time evolution. However, 
the extension of a natural vacuum free from IR divergence to 
the higher order in slow-roll expansion does not seem trivial at 
all. In this brief report we have also neglected the transverse- 
traceless metric perturbation, which can participate in the IR 
divergence. These issues will be discussed in our future pub- 
lication [7]. 
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